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THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS: 

1. A method for estimating the frequency of a single frequency complex 
exponential tone in additive Gaussian noise, comprising the steps of: 

5 performing the fast Fourier transform (FFT) on the tone; 

estimating the frequency as the frequency corresponding to the largest 
FFT output coefficient magnitude; 

computing a discriminant which is proportional to the frequency error in 
the initial frequency estimate using modified coefficients of the discrete Fourier 
10 transform (DFT) with center frequencies plus one half and minus one half of the 
FFT bin spacing relative to the initial frequency estimate; 

mapping the value of the discriminant into the estimate of the frequency 
error in the initial frequency estimate using a mathematically derived function; 
adding the estimate of the frequency error to the initial frequency 
15 estimate to get a first interpolated frequency estimate; 

computing a further discriminant which is proportional to the frequency 
error in the first interpolated frequency estimate using modified coefficients of 
the discrete Fourier transform (DFT) with center frequencies plus one half and 
minus one half of the FFT bin spacing relative to the first interpolated frequency 
20 estimate; 

mapping the value of the further discriminant into the estimate of the 
frequency error in the first interpolated frequency estimate using the 
mathematically derived function; and 

adding the estimate of the frequency error in the first interpolated 
25 frequency estimate to the first interpolated frequency estimate to get a second 
interpolated frequency estimate. 

2. The method according to claim 1 , wherein the first interpolated frequency 
estimate is in a region of relatively low noise induced frequency error. 

3. The method according to claim 1 or 2, wherein the method is 
30 implemented in computer hardware and/or computer software. 

4. The method according to any one of the preceding claims, wherein the 
method is utilised in communications, signal processing and biomedical 
applications. 

5. The method according to any one of the preceding claims, further 
35 comprising the steps of: 

iteratively deriving an interpolated frequency estimate, and 



i i 

WO 2004/005945 



1 I 

PCT/AU2003/000862 



31 

using the frequency discriminant, to obtain a more precise frequency 
estimate. 

6. The method according to claim 5, wherein the steps of iteratively deriving 
an interpolated frequency estimate and using the frequency discriminant are 

5 repeated until a fixed point solution occurs, where at this fixed point, the 
discriminant function has zero value. 

7. The method according to any one of the preceding claims, wherein the 
frequency discriminant is computed by: 

1 \fi\+\*\ 

k 

10 where, s = fr„ — 

N 

£ = fT s -:^L, anc j 
N 

p and a are the modified DFT coefficients defined by, 

1 N-l i2 - (kmm+ 2 ) 

^ = Y(k max+ i) = 2'( n )e N 

; 

15 and, a = Y0W -±) = £r(n)e ^"HT 2 - 

thus, the initial frequency estimate using the FFT, f 0 T s = ^s-and e « 0. 

N 

8. The method according to any one of claims 1 to 6, wherein the frequency 
discriminant is computed by: 

Y\p\ r +\a\ r 

20 where 0 and a are the modified DFT coefficients defined by, 

z n=0 

N1 (k—4> 

and,a = Y(k mflX -i) = 2r(fl)e J2an N 

*> n=0 

9. The method according to any one of claims 1 to 6, wherein the 
25 discriminant of frequency estimation error is computed by: 

P- Hll 2 -!*! 2 
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where 0 and a are the modified DFT coefficients defined by, 
P = Y(k raax + h = £r(n)e J2,m - r ~ 

1 n., < k ~4> 

and, a = Y(k max -i) = £r(n)e" J " 

Z n=0 

10. The method according to any one of claims 1 to 6, wherein the frequency 
discriminant is computed by: 

D-1**£^L] 
/3 + a 1 

where Re[.] is the real part and * denotes the complex conjugate, and p and a 
are the modified DFT coefficients defined by, 



10 



20 



^ = Y(k_ + I) = £r(n)e- j2,m » 



n«0 



and, « = Y(k max -I) = 2^)e* j2, ' a " * 



n— 0 



11. The method according to any one of preceding claims, wherein the 
frequency discriminant is computed by using more than two DFT coefficients. 

15 12. The method according to claim 11, wherein 2M+2 coefficients are used, 

N 

where 0<>M<, — -1 and the FFT coefficients are used in the frequency 

discriminant with optimal weighting coefficients obtained by using the concept 
of matched filtering is, 

E c _ i ,{Y[(k max +^+m)modN]-Y*[(k max -i-m)modN]} 
D =Re[ °°° ? I 3 

S {C ^V-] {YC(kmax + i +m ) modN l +Y '[( k m ax -^-m)modN]} 



where, 0<.M £ — -l, mod N indicates modulo N, 

2 

and, where, * denotes complex conjugate. 
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and, Y(k oax +-+m)andY(k ma3t -™-m) are the modified DFT coefficients given 
by. 

Y(k max +± + m) = £r(n)e » 

n=0 

5 

and, Y(k max -i-m) = 2r(n)e J 

■* n=0 

13. The method according to claim 11 or 12, wherein the frequency 
discriminant using more than two DFT coefficients us used in the last iteration 

10 to obtain additional frequency accuracy. 

14. The method according to any one of claims 11 to 13, wherein the 
frequency discriminant is computed by using more than two DFT coefficients 
and less or equal to all N FFT coefficients. 

15. The method according to any one of the preceding claims, wherein 
15 additional frequency accuracy is obtained by computing the frequency 

discriminant recursively until convergence for the frequency estimate is 
reached. 

16. The method according to claim 15, wherein convergence for the 
frequency estimate is reached after zero to three iterations, the number of 

20 iterations being dependent on the specific discriminant used and the signal to 
noise ratio. 

17. The method according to claim 15 or 16, wherein in any iteration, the 
frequency discriminant is computed using any one of the functional forms: 

AT.fr) ° 1 tan- 1 [}44^tan(^)] T 8 , Or 
25 AT m fr)=— j— [j ft 1 [" "l" 1 "!?] f s , Where y m isaconstant,r m > 0 , or 

Ar m fr)=-L[ '^H gm j 3T s , for r = i,or 

AFfrW— r ^ w t ~* < * m C ] for y=2 . 
m 4N L | J ff m | 2 +|a m p Jl8 ' 7 

1 8. The method according to claim 1 7, wherein y varies on each iteration. 

19. The method according to claim 15 or 16, wherein in any iteration, the 
30 frequency discriminant is computed using: 

Af m (r)=-l-Re[ ^" g ? ]f.. 

2N Pm + a m 
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where, Re[J denotes the real part and * denotes the complex conjugate. 

20. The method according to any one of claims 17 to 19, wherein the 
frequency incremental shift, Af m (r), is related to the previously defined 
frequency discriminant, D, by, 

21. The method according to any one of claims 15 to 20, wherein the 
frequency discriminant is driven to zero input and output values by either 
modifying the frequency of the DFT coefficients or frequency translating, the 
signal. 

10 22. The method according to any one of claims 15 to 21, wherein signal 
frequency translation is achieved by multiplication of the signal by a locally 
generated complex exponential signal. 

23. The method according to claim 22, wherein frequency multiplication of 
the signal is implemented with a standard hardware, software, or combination 

15 hardware/software FFT. 

24. The method according to claim 23, wherein the hardware/software FFT 
is highly optimized for at least one processor operating as a system. 

25. The method according to any one of claims 15 to 24, further comprising 
the step of scaling the frequency estimate during recursion, to save multiplies. 

20 26. The method according to claim 25, further comprising a final step of 
multiplying the scaled frequency estimate f m+1 T 8 with the sampling frequency f 8 
to remove the scaling from the frequency estimate. 

27. A frequency estimation software program for estimating the frequency of 
a single frequency complex exponential tone in additive Gaussian noise, 

25 wherein the frequency estimation program has functionality to perform the 
method according to any one of claims 1 to 26. 

28. A computer system programmed to perform the method according to any 
one of claims 1 to 26. 

29. The computer system according to claim 28, wherein the hardware 
30 includes a DSP processor chip. 



